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Abstract: We perform a general study of the primordial scalar non-Gaussianities in multi- field 
inflationary models in Einstein gravity. We consider models governed by a Lagrangian which is a 



lO \ general function of the scalar fields and their first-order spacetime derivatives. We use SN formalism to 



relate scalar fields and the curvature perturbations. We calculate the explicit cubic-order perturbative 
action and the three-point function of curvature perturbation evaluated at the horizon-crossing. Under 
reasonable assumptions, in the small slow-varying parameters limit and with a sound speed c s close to 
one, we find that the non-Gaussianity is completely determined by these slow-varying parameters and 
some other parameters which are determined by the structure of the inflationary model. Our work 
generalizes previous results, and implies the possibility of the existence of large non-Gaussianity in 
model constructing, and it would be also useful to study the non-Gaussianity in multi-field inflationary 



models which will be constructed in the future. 
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1. Introduction 

Inflation has been a very successful paradigm for understanding the evolution of the very early universe 
[PJ. It not only naturally provides a way to solve flatness and horizon problems, but also generates 
density perturbations as seeds for the large-scale structure in the universe. Inflation is most commonly 
discussed in terms of a potential energy which is a function of a single, slowly rolling scalar field. Such 
models generically predict an almost scale-invariant spectrum and an almost Gaussian distribution of 
adiabatic density perturbations on super-Hubble scales [§]. These generic predictions are consistent 
with recent cosmological observations of the temperature anisotropy and polarization of the Cosmic 
Microwave Background ||. 

The measurement of primordial perturbations provides increasingly precise determination of the 
spectrum index n s and tensor-to-scalar ratio r, which come from the power spectrum, in other words, 
the two-point correlation function of the curvature perturbation. However, there are many alternative 
models of inflation are able to give the similar predictions, and thus there is still considerable ambiguity 
in constructing the real inflationary model. 

In contrast, the non-Gaussian component of the scalar perturbations is characterized by the cor- 
relation functions beyond two-point, e.g. the three-point function of the fluctuations, which is a 
nontrivial function of three variables and will provide us more information beyond the power spec- 
trum. Furthermore, as described above, the non-Gaussianity of distribution of primordial fluctuations 
predicted by the simplest model of inflation is well below the current limit of measurement ffl, 0, pi ||] . 
Therefore, any detection of large non-Gaussianity would be a significant challenge to our current 
understanding of the early universe. 



Indeed, there is the possibility of the presence of large non-Gaussianity with /nl S> 1 H §• The 
latest observational bound on the three-point function of the primordial curvature perturbation £ from 
the WMAP 5-year data states that the local non-gaussianity parameter /nl is limited to the values 
—9 < /nl < m Hi- I n ^ ne next few years, with improved experiments like the Planck satellite, we 
will measure the CMB anisotropics to an incredible resolution at |/nl| < 5. 

On the theoretical side, the level of non-Gaussianities can be generally calculated analytically in 
typical single or multiple- field inflationary models (see ||] for a review). The three-point function for 
the standard single-field inflation models with a canonical kinetic term was performed in an elegant 
and gauge invariant way in [y, 0, and the corresponding /nl is O(10~ 2 ) and too small to be detected. 
In order to find possibly detectable large non-Gaussianities in single-field cases, models with non- 
standard actions have also been studied [21, 22, 2£, 27, fZlJ] . 

In the past few years, intensive efforts have been devoted to connecting string theory and inflation 
pQ , O , |12] , 13, |lj]. Indeed, in models descending from the low-energy limit of superstring theory or 
string compactification, there are many light moduli fields describing the higher-dimensional degrees of 
freedom, which could play a role during inflation. Therefore, it is interesting and of physical significance 
to generalize the above analysis of non-Gaussianities in single-field models to the multiple-field inflation 
models M, M, M, p], |1 P, M, M EJ7L pL |T ' 



gg, gg, gg g§ g|, g§ g§ g7|]. Most of the previous 
works focus on the models with standard canonical kinetic terms for the scalar fields involved. It is 
shown that in models satisfying the slow-roll conditions, /nl is also of order slow-roll parameters as in 
the standard single-field case. Models with non-canonical kinetic terms have also been investigated, 
in particular the effective multiple- field DBI- inflation [fig, g6|, fiT]]. Models with non-trivial field space 
metric Gjj have also been considered flcfl . 

In these inflationary models inspired by string theory or other theory beyond the Standard Model, 
there are several significant features different from those in the standard single field models, where 
the Lagrangian is of the form C = \f—g [— \{d(j)) 2 — V((f))] = ^—g(X — V). Firstly, the scalar field 
potential is not necessarily the only degree of freedom in model building, in contrast to that in the 
standard single field models. Indeed, in these string inspired models, it is generally expected that 
deviation from the standard kinetic term X ~ {d(j)) 2 of the scalar field action would arise ]l~5| . One 
would generally expect loop corrections in the quantum theory to generate operators in the Lagrangian 
that are proportional to higher-order derivatives X 2 , X 4 , and so on. Indeed, if the energy scale of 
renormalization is of order Planck scale M p \, such higher-order interactions would be suppressed and 
be negligible, and thus the canonical kinetic term would give a good approximation. On the other 



hand, if the inflation scale is close to the GUT scale, these ultraviolet corrections might be significant 
and of considerable relevance [16]. Such models with non-canonical kinetic terms have been considered 
previously by a number of authors [16, 17, [l8|, [l^, |20[| . For our purpose, the deviation from a canonical 
Lagrangian would be an important source of possible large non-Gaussianities. 

Secondly, the presence of many scalar fields during inflation would affect the generation of pri- 
mordial perturbations. In single field inflation, it is well-known that the curvature perturbation is 
conserved on super-horizon scales. However, in multi-field inflation models, the presence of multiple 
light fields will lead to the generation of non-adiabatic field perturbations during inflation. The curva- 
ture perturbation in multi-field inflation can generally evolve after Hubble exiting, due to the presence 
of entropy perturbations, which can be 'sources' of curvature perturbation. This non-trivial evolution 
of the overall curvature perturbation would be another source of detectable non-Gaussianity. A general 
description of calculating the three-point functions of the fluctuations in multi-field models was pre- 
sented recently in [35], which emphasized this super-horizon evolution of curvature perturbation. It is 
also argued in |j(| that in some cases, for example in the curvaton scenario |^, 55, ^, 57, 5^, [59, 72, |6l|l , 
there is the possibility that the large-scale superhorizon-effected components would dominate the pri- 



mordial non-Gaussianity, rather than the components from the microscopic fluctuations. However, 
in the absence of any estimation for the microphysical contributions to the three-point functions, the 
above analysis is incomplete, and thus a precise investigation of non-Gaussianities generated from the 
fluctuations at the horizon-crossing is needed. 

In this work, we consider a very large class of multiple-field inflation models, in which the dynamics 
of scalar fields is described by a Lagrangian P(X, (j) 1 ), where P is an arbitrary function of J\f scalar 
fields and their kinetic term X = —^(Vfj} 1 ) 2 |25[ |. This form of action includes the standard choice 
P = X — V as a special case, and can be viewed as a generalization of the Lagrangian of k-inflation 
Ijlql to the cases of multiple scalar fields. Low-energy limit of string theory may also lead to such 



an action, such as Dirac-Born-Infeld action 23, 24], and its multi-field extension studied recently in 



Ea , 46]. In general, one may consider theories including arbitrary number of higher-order derivatives 
in the action, such as P(4>, d(f>, d 2 <p, ■ ■ • ). However, if the energy scale of inflation process is much 
lower than M p \, the contributions from these higher-order derivatives will be suppressed and thus can 
be neglected. Therefore, we consider models where P contains only X. Due to the same reason, we 
consider Einstein gravity. 

We assume that these scalar fields generate the density perturbations. We expand the general 
multi-field Lagrangian to cubic order of the perturbations 8(f) 1 , and use SN formalism (48|, |49], p^H, 



|5l| ] to relate curvature perturbation £ and these multiple scalar perturbations 8<fr . We calculate 
the scalar three-point function of the curvature perturbations following Maldacena et. al. M. To 
control the calculation, we define some small slow-varying parameters, which can be viewed as the 
generalization of the standard slow-roll parameters. In our formalism, although the Lagrangian is very 
general, we assume the effective speed of sound c s is almost close to one. In the limit of small slow- 
varying parameters, under some reasonable assumptions, we finally find that the non-Gaussianity in 
these models is completely characterized by the slow-varying parameters and some other parameters 
determined by the concrete structure of P(X,(j) 1 ). In particular, our result shows the possibility of 
the presence of large non-Gaussianities, due to the deviation from the canonical Lagrangian. 

The paper is organized as follows. In the next section, we setup the general multi-field inflationary 
models, and derive the background equations of motion. In order to control the calculation and 
analysis the solution, we define some slow-varying parameters in this general context. Then we give 
a brief review of calculating non-Gaussianities during multi-field inflation based on 5N formalism. In 
section 3, we develop the second-order theory for the linear perturbations. We estimate the power 
spectrum in the limit of small slow-varying parameters. In section 4, we calculate the exact cubic 
perturbation action for this most general multi-field inflationary model. In section 5, we perform a 
general calculation of the three-point function, which represents the central result of this work. Finally, 
we make a conclusion in section 6. 

We work in natural units, where c = h = M p \ = 1, and M p \ = (8-7rG) -1 ' 2 is the reduced Planck 
mass. 

2. Setup 

2.1 Background 

We consider a large class of multi-field inflation models, which are constructed from a generic set of 
M scalar fields {(j) 1 ,1 = 1,2, •• • ,J\f} coupled to Einstein gravity. The action takes the form: 



S= j d 4 x^j ^R + P(X,^ 



(2.1) 



with kinetic term 

X = —Gug^d^dutf , (2.2) 



2 



where g^ v is the spacetime metric with signature (— , + + +) and Gjj = Gu(4>) is the AA-dimensional 
field space metric, and P is an arbitrary function of X and </>'s. 

The energy-momentum tensor of the scalar fields takes the form 

Tl w = Pg ^u + Px G IJ d^(t) I d u ( t ) J , (2.3) 

where Px denotes the partial derivative of P with respect to X. In order to consider the back- 
ground (unperturbed) dynamics, we suppose that the universe is homogeneous, with a flat Friedmann- 
Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (t)dx l dx i , (2.4) 

where a(t) is the scale factor and H = a/a is the Hubble parameter. Under this assumption, the 
energy-momentum tensor of the scalar fields reduces to that of a perfect fluid, with energy density 

p = 2XP X - P , (2.5) 

and pressure P = P(X, 4* 1 ). 

The equations of motion for the scalar fields derived from ( |2.1| ) are [25] 




, (2.6) 

where Pj denotes the derivative of P with respect to (p 1 , and Fj K is the Christoffel symbols associated 
with the field space metric Gjj . The equations of motion of the gravitational dynamics are Friedmann 
equation 

T 2 P _ * , 



H Z = ^-(2XP X -P), (2.7) 



and the continuity equation 



p = -3H(p + P) = -QHXPx ■ (2.8) 

The combination of the above two equations gives another useful equation 

H = -XP X ■ (2.9) 



It is convenient to define an effective speed of sound c s [16|, as 



P,x Px + 2a Pxx 

Note that in models with canonical kinetic terms P = X — V, c s = 1. 

In general one may consider models with an arbitrary field space metric Gu^ 1 ). In this work, 
we focus on the case where Gu = 5u and thus Tj^ = 0, i.e. the field space is flat. This choice of 
action already covers a large class of multi-field inflationary models with non-canonical kinetic terms. 
The linear perturbations of models with an arbitrary metric Gjj have been investigated in [pBj ], 

2.2 Slow-varying parameters 

For general function P(X, </> ), it is difficult to solve the equations of motion for the scalars (|2] 



analytically. In order to capture the main physical picture and investigate the evolution of the system, 
the idea is to define some small parameters to control the dynamics, and to find solutions perturbatively 
in power expansions of these small parameters. 



In standard single field slow-roll inflation, this condition is achieved by assuming that the inflaton 
(j) is rolling slowly in comparison with the expansion rate \(j)\ <C H. Similarly, in general multi-field 
inflation models, it will prove convenient to introduce a dimensionless "slow-varying" matrix as 



2H 2 

where 



£ » = £i|V v, (111) 



When there is only one inflaton field involved, e IJ reduces to the standard single field slow-roll pa- 
rameter e = —H/H 2 , which can also be expressed as 

e = tre /J = G u e IJ = --^. (2.13) 

It proves useful to decompose e into two new parameters ex and e^ [26], which measure how the 
Hubble parameter H varies with the kinetic and potential parts of scalar fields (j) 1 respectively 

e = e x + H = -^-X - ^ . (2.14) 

As in the single field models where we may define rj = jj-, here we define another slow- varying 

i.i 



matrix r/ as 



JJ e 



V = jji , (2-15) 



which can be written explicitly 



i j = 2e u _ paVV + V^ + pxVV (2 16) 

2HH 

The matrix rj IJ generalizes the slow-roll parameter 7] in single-field inflation models. Note that with 
a flat target space metric Gu, we have r\ = trr] IJ = e/eH as expected. 

Due to the generality of the fuction P(X, (p 1 ) and the complexity of the scalar equations of motion 
( |2.6D , the relations between Pj and the rolling of scalar fields (j) and <p are complicated. And thus 
we may introduce another set of parameters defined as 

h = - vfe ■ (217> 

which can be viewed as the analogue of the standard slow-roll parameter of the form e = \{V' /V) 2 . 
In models with standard canonical kinetic terms, it is easy to show that \e \ ~ \e \ as expected, and 
thus two slow-roll parameters e and rj are enough to control the theory. However, for general function 
P, there is no simple relation between these two parameters. In order to proceed, we expect that both 
e and e are of the same order, and assume that P t x/{HP,x) can be negligible. Therefore, from the 
scalar equations of motion ( |2.6| ), it is easy to see that e ~ —e as expected. Similarly another matrix 
is defined as 

this is the analogue of rj = V" fV ~ V" /3H 2 in the single field case. 
Furthermore, we define the dimensionless parameters 



1, « = ^7, (2.19) 

s 



c 2 c,H 



where u measures the deviation of the effective sound speed c s from unity, and s measures the change 
speed of c s . In models with canonical kinetic terms, u = s = 0. 

These parameters generalize the usual slow-roll parameters, and in general depend on the kinetic 
terms as well as the potential terms. For generic theories we expect that |e|, \u\, \s\ <C 1 and 1 



e 1J , V 1J ,fj 1J ~<D(e/N) , (2.20) 



and thus e 1 ~ e ~ I \fejM I , here M is the number of the scalar fields. Furthermore, for those 
models with non-vanishing "cross" derivatives, i.e. P t xi, P,xxi,P,xu / 0, we assume that the 'X'- 
derivatives of these slow- varying parameters also satisfy some smallness relations 

~i P,X ~j ~i P,x ~t ~tj P.x ~tj , n „-, x 

> x ~ ip ' ' xx ~ iF ' v ' x ~ jp v ■ ( ^ 

In general, the validity of these conditions depends on the explicit forms of the models. In this work, 
we do not try to find explicit models which satisfy these slow-varying conditions. It is also addressed 
in [p6| , ^7]] and etc., in the presence of a non-canonical kinetic term, the smallness of these slow- varying 
parameters does not imply that the inflation itself is slow-rolling. 

In order to simplify some of the derivations that follows, it will prove useful to define two param- 
eters which are combinations of derivatives of P with respect to the kinetic term X 2 

S = XP.x + 2X 2 P. XX = ^r , (2.22) 

A = X 2 P tX x + ~X 3 P X xx ■ (2.23) 

Here S is of order 0(e), but for general P(X, <fr ) where Pxi ^ 0, such as K-inflation or DBI-inflation, 
there is no simple relation between A and the above slow- varying parameters. 

2.3 5N Formalism 

In this section we make a brief review of 5N formalism [48, 49, 50, 51|, which is proved to be a powerful 
technique to calculate the curvature perturbation in a variety of inflation models, especially in the 
multi-field models. 

The idea of 5N formalism is to identify primordial curvature perturbation £ with the perturbation 
of the local expansion. Starting from a flat slice at some initial time t- l} the local expansion iV(i,t;,x) 
at some final time t is defined as 

N(t,ti,x) = f dt'H(t',x) , (2.24) 

Jti 

where H(t,x) is the local Hubble expansion rate due to the perturbations. Then the primordial 
curvature perturbation can be expressed as 

C{t, x) = N(t, t u x) - N Q (t, ti) = SN , (2.25) 



1 In this work we consider the case with c s is very close to one by assuming u ~ 0(e). In general, one may consider 
models with an arbitrary c s , however, as addressed in [149, M where the multiple DBI inflation is investigated, the 
multi-field effect is suppressed in the limit of c 3 <C 1, which may be a possible source of large non-Gaussianity in single 
field models. We will leave this for a future investigation. 



From (P 1 1C|),(2.22 )-( 2.23), we can extract some useful relations for later convenience 



X 2 P, XX = E(l - <2)/2 , X s P,xxx = 3A/2 - 3S(1 - c 2 s )/4 , \-l = 2X?' X A 



cl P 



X 



where No(t, U) is the background (unperturbed) expansion, which is related to the background Hubble 
expansion rate Ho(t) as 

No(t,ti)= j dt'H (t r ). (2.26) 

If we take t\ at the time of horizon-crossing during inflation denoted by £*, then N(t,t*,x) becomes a 
function of the scalar fields evaluated at horizon-crossing. Then £ can be expanded as 

C(t,x) = ^iV i/ (£)^(x) + i^7V, 7J (t) ( 5^(x) ( 5^(x) + ... , (2.27) 

/ I, J 

Note that £(£) is just the primordial adiabatic curvature perturbation if we choose t well after the 
reheating process. 

After going to the momentum space, we have 

C(k) = Nj5^(k) + ^NuW 1 * <^ J ](k) + • • • . (2.28) 

The two-point and three-point functions of C, can be expressed in terms of the two and three-point 
functions of the scalar fields fluctuations 5<p as 

(C(ki)c(k 2 )> = Jv;zJv; 1 7<<ty / (ki)<ty J (k 2 )> + • • • , (2.29) 

(C(ki)C(k 2 )C(k 3 )) = iV,/iV,jiV^(^ / (k 1 )(5 ( / ) J (k 2 )^(k3)) 

+ iiV i/ iV i j^ L (^ / (k 1 ) ( 5 < /» J (k 2 )[^*^ L ](k3)) + perms+--- , (2.30) 



where * denotes the convolution product. In Section ^^ we can see that the two-point functions for 
the scalar fields satisfy 

2?r 2 
(^ 7 (k 1 )<50 J (k 2 )) = (27rf5 2 (k 1 + k 2 )G /J -^A2 , (2.31) 

and in Section [| we will show that the scalar three-point functions can be written in the form 

(^ / (k 1 )^ J (k 2 )^(k 3 )) = (2^) 2 <5 3 (k 1 +k 2 + k 3 )-^-\Al\ 2 A IJK , (2.32) 

where A 2 is the power spectrum of a massless scalar field in de Sitter space. The principal result of 
this work is the momentum-dependent function A. (ki,k2,ks) given in ( 5.15| ), which contains the 
information of the amplitude and shape of the non-Gaussianity 

In order to connect the above analysis with the observations, the non-Gaussianity measured by 
the three-point functions must be expressed in terms of an experimentally relevant parameter. A 
common choice is the non-linearity parameter /nl defined as 

C = Cg + jj/NLC g 2 , (2-33) 

which denotes the departure of £ from a Gaussian random variable ( g . The power spectrum and 
bispectrum of £ are defined in terms of the two and three-point functions respectively as 

(C(ki)C(k 2 )) = (2^) 3 5 3 (k! + k 2 )P c (fcx) , 

(C(k 1 )C(k 2 )C(k 3 )) = (27r) 3 <5 3 (k 1 + k 2 + k 3 )i? c (A;i,A ;2 ,fc3) , 

then B^ is related with P^ in terms of /nl as 

B C = pNL(h,k 2 , fc 3 )[P c (A: 1 )P c (A ;2 ) + P ( (k 2 )P c (k 3 ) + P c (A&)P c (*i)] . (2.35) 



To relate /nl with the three-point function of £, we write 

(C(k 1 )C(k 2 )C(k 3 )) = (2vr) 3 5 3 (k 1 + k 2 + k 3 ) jpg \Al\ 2 Ac , (2.36) 

then /nl can be written as 

/NL= ^E^r (2,37) 



From ( p. 29 ) and (2.30), A( can also be written as 

„4 C = N I N J N K A IJK + G IK G JL N I N J N KL ^ kf . (2.38) 

i 

From the above relations, the non-linearity parameter /nl can be expressed in terms of N j and the 
momentum-dependent function A , up to leading orders, as 

5 NjNjN k A ijk b G IK G JL N yI N yJ N yKL 
/NL 6 (G IJ NjN :J y £. fcf + 6 (G^iV/A j) 2 + " ' ' [ ' 



where '• • • ' denotes the remaining cross terms from ( 2.30J ) which we have neglected together with other 



higher-order terms. This expression was first derived in paL 

Furthermore, in order to calculate the primordial power spectrum and the non-linearity parameter, 
we need to know the derivatives of the number of e- folding N with respect to the scalar fields N j, 
Njj, etc. It is easy to show that dN = —dlnH/e, and thus 



-V/"-\/^7 + - • (2-40) 



3. Linear Perturbations 



3.1 ADM formalism and the constraint equations 

In single field inflation models, we have two different gauge choices. One is the gauge where we 
consider the curvature scalar on uniform density hypersurfaces, defined by 5(p = 0; the other is the 
spatially flat gauge, where we choose the uniform curvature slicing, and the spatial part of the metric 
is unperturbed. The physical degrees of freedom of the perturbations are completely described by the 
perturbations of the metric in the first gauge, and only by the perturbations of the scalar fields in the 
spatially flat gauge. In single-field case, both two gauge choices are natural. However, in the case of 
multi-field models, the first gauge is no longer possible, and the only natural choice is the spatially 
flat gauge, where the physical degrees of freedom are perturbations of the scalar fields. 
It is very convenient to work in the ADM metric formalism 



ds 2 



-N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) , (3.1) 



where iV is the lapse function and N l is the shift vector. The ADM formalism is convenient because 
the equations of motion for N and iV* are exactly the energy and momentum constraints which are 
quite easy to solve. 

Under the ADM metric ansatz, the action becomes 3 

S = I dtd 3 xVhN i?— + p\ + j dtd z x^- {EijE ij - E 2 ) , (3.2) 



3 Here and in what follows, the spatial indices i, j are raised and lowered using hij. 



where h = det/ij,- and the symmetric tensor E^ is defined as 



E i j = -[h i j-V i Nj-VjN i ) , (3.3) 



and E = tiEij = h 13 Eij. R^> is the three-dimensional Ricci curvature which is computed from the 



metric hij. The kinetic term X now can be written as 



X = 2iv^ G/j7rV ~ %W<%0 J > ( 3 - 4 ) 

with 

vr 7 = ft - N i d l cp I . (3.5) 

The equations of motion for N and N l give the energy constrain and momentum constraint 
respectively 

P ~W2 i E » EiJ ~ E2 + ZPxGutt 1 ^) = , (3.6) 

V, (i <E{ - E5{)} = ZZ-GuJOv? . (3.7) 

In spatially flat gauge, we have hij = a 2 (t)5ij and thus w 3 ' = 0. The unperturbed flat FRW 
background corresponds to N = 1, N l = 0. The scalar fields on the flat hypersurfaces can be 
decomposed into 

I (t,x) = 4(t) + Q I (t,x), (3.8) 

where (J)q are the spatially homogeneous background values, and Q are the linear perturbations. 
In what follows, we always drop the subscript '0' on <pQ and simply identify cj) as the unperturbed 
background fields. In order to study the scalar perturbations of the metric and the scalar fields, we 
may expand N and N % as 

N = 1 + ai + a 2 + ■ ■ ■ , , N 

(3.9) 

jv* = d l p = d l (/?i + & + ■ ■ ■ ) , 

where a n ,(3 n are of order 0{Q n ). One can plug the above power expansions into the constrain 
equations of N and N l ( |3,6| )-( |3"77| ) to determine a n and (3 n . To the first-order of Q the solutions are 



«i = ^iQ 1 . ( 3 - 10 ) 



and 



** = W 



■^fhQ 1 - 2XP XI Q J + PjQ 1 + ^ (^- - 3H 2 \ ^Q 1 



(3.11) 



where d 2 = 5 lJ didj. Fortunately, it turns out that in order to expand the effective action to order 
0(Q 5 ), in the ADM formalism we do not need to compute N and N l to order 0(Q 3 ), since they 
must be multiplied by dL/dN or dL/dN l which vanish due to the constraint equations. Also in the 
present case, terms of order 0(Q 2 ) in N and N l drop out of the third-order effective action, and thus 
( |3.10| )- (|3.11| ) are sufficient for our purpose. Furthermore, it is easy to see that ct\ and f3\ are both of 
order C(e 7 ). 



3.2 The second-order theory 

In the Appendix [A], the general form of the expansion of the action to the cubic-order of Q has been 
developed. From flA.ll ), the second-order action can be written as 4 [25] 

S {2) = \j dtd 3 xa 3 [(P,xGu + Pxx4>ih) Q'Q J + 24> I P XJ Q I Q J - P^Gu&Q 1 ' d % Q J - MuQ'Q 3 

(3.12) 
where Mij is the effective mass-matrix 

2XP X : „ XP3 X (. l\ : : Id 



2XP X ■ Xp x ( 1 \ • • Id 

Mu = -Pu + -^<f,,P XJ + -^{i-^ «, - -,- 



^ 1 + 5 IW ' 



(3.13) 



( |3.12| ) is exact and is valid for arbitrary scalar fields dynamics. 

In general, one may decompose the multi-field perturbations into an adiabatic mode and N — 1 
entropy modes, since the background trajectory specifies a special inflaton direction. Here in this 
work, under the assumption that c s is close to unity, we may take a different method, to analyze all 
these perturbations in a unified and more simpler formalism. To solve the second-order action ( 3.12J ), 
we define a symmetric matrix ajj which satisfies 

aiKORj = a 2 P tX Au , (3.14) 

with 

Au = Gu + ^^M'j = Gu + -e u , (3.15) 

a> P,x e 

where a prime denotes the derivative with respect to the conformal time r\ defined by drj = dt/a. 
Here Au denotes the deviation from the standard canonical action, which reduces to Au = Gu 
when considering the standard kinetic terms, where u = 0. Then we define the new Mukhanov-type 
variables 



ii 



a IJ Q J , Q 1 = ( a - 1 ) IJ u J , (3.16) 



or in matrix form u = aQ and Q = a~ 1 u. This treatment is an analogue of that in the single-field 
case. For standard single-field inflation models with a canonical kinetic term, au reduces to a, and 
uj = ajjQj reduces to the usual rescale relation u = aS(j), and thus the matrix ajj can be viewed as 
the generalization of the parameter a in single-field case, where we may rescale the field as u = aS(j>. 
In terms of these new variables u , after changing into conformal time 77, the second-order action 
( |3.12| ) can be rewritten in a matrix form as 

£(2) = f dr]d 3 x \u' T u' + u' T i^a'a- 1 + 2a 2 a~ 1 Ba~ 1 ) u - a 2 A- 1 d i u T d i u 

J ~ (3.17) 



+u (—2a a a'a Ba —a a Ma +a a'a'a ' ) u\ 



a a 



l„,-i 



where we have defined Bjj = 4>'jP^x.j for short. In deriving fl3.17j ) we used the fact that 
A~ 1 /a 2 P j x- The equations of motion for the scalars can be derived from fl3.17j ) and written in a 
compact form 

u" + la' 1 a' - a'a' 1 + a 2 a - l (B - B T )a- 1 } v! 

f3 18") 
+ [A' 1 k 2 + a 2 a' 1 (a 2 M + 2HB + B' - (B - B T )cT V) a' 1 - a" a' 1 ] u = 0, 

where TL = a' /a. Note that ( ^.18[) is exact and no approximation is made. For standard single field 
models, ( |3.18| ) reduces to the well-known result u" + (k 2 + a 2 m 2 — a" / a)u = as expected. From ( |3.15| ) 



This expression should be compared, for example, with eq. (44) of [B6| , to which it reduces in the multi-field models 
with canonical kinetic terms where P = X — V . 



it can be seen easily that in the single- field case, the matrix Au reduces to a pure number A = 1/c^, 
and thus the first term in the second line, which is proportional to the wave number k, reduces to 
c 2 k 2 as expected. 

We would like to make some comments here. Prom ( [3.15 ) we can see that the unperturbed back- 



ground fields velocities 4>'j, represent a special direction in field space when considering perturbations 
of the fields. This is just the so-called adiabatic direction which has been introduced in [52] for multi- 



field inflation models with canonical kinetic terms. The decomposition into adiabatic and entropy 
modes is equivalent to a 'local' rotation in the space of field perturbations, 

Q n = e^Q 1 , Q 1 = e{Q n , (3.19) 

the rotation matrix e" is just the projection of a new set of basis {e n } on {<pi}- The first vector is 
specified as 

e\ = -p= . ( 3 - 2 °) 

V2X 

which is just the unit local adiabatic vector. Note that the rotation e" is defined locally, and dependents 
on the background trajectory. The new field Q 1 is just the adiabatic mode while other modes Q n (n = 
2, • • • , jV) are entropy modes. After the local rotation, i.e the decomposition into adiabatic and entropy 
modes, it is easy to show that Au = Gjj + (l/(? s — l) e}ej, and 

AuQ'QJ =r -L (Q 1 + ZlQP) (q 1 + Zffi) + Y, (& n + ZpQ p ) (Q m + ^ m Q 9 ) , (3.21) 



where we denote Z™ = efe^ for short. This decomposition clearly shows that the adiabatic component 
of the perturbations Q 1 obeys a wave equation where the propagation speed is the sound speed c s , 
while the entropy modes of the perturbations Q n (n = 2, • • • ,J\f) propagate with the speed of light 



c = 1. This property was first pointed out in the studying of two-field DBI-inflation [45, |4q| , but here 
we see that it turns out to be a generic feature for general multiple field inflation models. 

The decomposition into adiabatic and entropy modes is a powerful tool to analyse the equations 



( |3.18| ). In general, however, (3.18) is a set of coupled equations and rather complicated to solve, even 



in the case of canonical kinetic term. The idea is to take slow-roll approximations as in the standard 
slow-roll inflation models. This is the standard approximation in estimating the amplitude of the 
perturbation power spectra, i.e. the two-point functions. For our purpose, it also supplies a unified 
treatment with all these perturbations Q , rather than decomposing them into adiabatic and entropy 
modes. 

In single-field models, the mass term can be expressed as a combination of slow-roll parameters, 
and thus be neglected to leading-order in slow-roll approximation. Similarly, the mass-matrix defined 
in fl3.13|) can also be written in terms of slow-varying parameters as 

jpM u = -3P jX ( 1 + -g J eu + 3f)u 

12 ( P ' x ~ ^~ \ P > x (i X \ ( p A ( lN lf ^ 2P >* 



= -6P x eu + 3fi I j + 0(e 2 ). 

(3.22) 
We also note that Bjj = (f>jP t xj ~ C( e )- Since we assume that c 2 s departs from unity by a quantity 
that is first-order in slow-roll, i.e. {l/c 2 s — 1) ~ 0(e), Au can be written as 

Au = Gu + 0(e) , (3.23) 



and thus (.4 l )u = Gjj — 0(e). Furthermore, from ( |3, 14 ) it is easy to show that the last term in 
(|3.18|) can be written as 

a "cT 1 = — (l + 0(e)). (3.24) 

a 

Bring all these considerations together, the equations of motion to the lowest-order of slow- varying 
parameters reduce to a very simple form which are just M decoupled de Sitter-Mukhanov equations 

u i" + ( k 2_?pj u i = 0: (3.25) 

where a" /a = 2/rj 2 + 0(e). 

The solutions to ( |3.25D are standard, we then find the Q 1 two-point functions as 

(Q / (r ?1 ,x 1 )Q J (r ?2 ,x 2 )) = G /J G*(7 7 i,xi;r ? 2,x 2 ) , (3.26) 

while G* represents 

H 2 f (1 + ik m )(l - ik V2 )e-^-^ , m>m 

GJm,m,k) = — x< .,, , (3.27) 

M ) 2k3px | (l _ ,^ )(l + . k ^ )e+lk{m _ m) ^ m<r)2i 

where we have chosen boundary conditions so that G* behaves like flat space propagator at very early 
times, when the perturbation modes are deep inside the horizon. This corresponds to the Bunch-Davies 



vacuum [53|. The power spectra for the scalar fields on large scales can be read easily 

(Q^kOQ^ka)) = (2vr) 3 <5 3 (k 1 + k 2 )G /J ^- , (3.28) 

and therefore the dimensionless power spectra are 5 

A - = i^k ■ (3 - 29> 

4. Non-linear Perturbations 

4.1 General form of the third-order action 

In this section we turn to the central calculation of this work, the third-order piece of the coupled action 
( |3.2[) . The general form of the expansion of the action to the cubic-order of scalar fields perturbations 
Q 1 is developed in Appendix [A] 

s® = I dtd'xa 3 (V 3 ) + ai pM - ai nw + ofnW - « 3 n(°)) , (4.1) 



5 This result coincides with (97) and (105) of J25J in the c s — » 1 limit as expected, while in |Eq| the adiabatic and 
entropy spectra are obtained respectively with an arbitrary c s . 



where a.\ is given by ( 3,1C| ) and pw ) n' n ' can be found in Appendix [A]. After a straightforward but 
rather tedious calculation, we find 



S( 3 ) = / dtd 3 xa 3 

( 3AP 3 , 



P 2 
,x 



3AP 3 X , . . l, . „ j 



IAJ/^K 



Jh^Gu^k ~ 4gs 2 c 4 ^ J ^ J ^ + 2 GljP,XK + 2^ I< ^ jP ' XXK j Q Q Q 



+ 






l ^ H2( p Vifai>K - ^MjPxk - ^hhP,xxK + \4>iP,xjk ) g 7 Q J Q* 



,A' 



P 3 , 



P 



4ff 2 S <t>I<t>jP,XK ~ ^p( 2A - S + 3H Z )(f)l(t>J(t)K + -^24>l4>jP,K + -^jrP,IJ(f>K 



P 2 

r X 



2„4 



+ 



E 2 c 



Sc? 



2jj2 9l<t>jP,XXK ~ -^4>lP,XJK + gP/JX ) Q^Q' 



\InJr\K 



v/^~>i f) r\J dK 



+ (^PxxGu^k + ^PxxxMAk) &Q j Q k - -^GufaffQ'W'Q' 
P 2 

- PxQi&PidiQ 1 + ■£i>ihQ I &p\diQ J - P,xxhhQ I d i (3 1 d i Q J 

Zri 



P,xxZc 2 



F ^LG IJ ^ K d i Q I d i Q J Q K + 

1 



nhQ I d l (3id l Q J + Pxx "'"" : 



irJn nJr>K 



2H 



G IJ 4>Kd'Q 1 d l Q J Q 



P 



P,xikjQ I 9 L (3id l Q J - -G I jP,x K d i Q I d i Q J Q K - ^iQ 1 (^WyA - (d%/3x) 2 ) 

' (4.2) 

No approximation of small slow- varying parameters has been made in deriving ( |4.2j ), and thus it is 
exact. 



4.2 Slow-varying limit 



In order to proceed, we restrict (4.2) to the leading-order of slow-varying parameters. This is because 
of not only the complexity of the full cubic-order action Q4.2] ), but also the observational constraints. 
Therefore the third-order action can be written in a much simpler form 6 



5 (3) 



dtd 3 xa 3 



SIAJAK 



\p,xxGij^K + \p,xxx»iO., IK j ()'()■>()' 
P 2 \ 

r ,x ■ * 



+ I \gijP,xk - ^§ g u^k) Q 1 Q j Q k - ^^Gij4> K d i Q I d % Q J Q K 



(4.3) 



X Gu4> K + 1 -G u P,xk J d i Q I d i Q J Q K - PxQid^idiQ 1 



/p2 

~ \7h~' ut " ' 2 

Pi is defined in ( p.ll ) and is also given to the leading-order by 

a 2 P 



d 2 (3i 



,x 



2H 



uQ 1 . 



(4.4) 



Now it proves most convenient to take integral by parts to eliminate the last term containing di(3\ 
in (fOj). It is easy to show that 



dtd^xa^PxQid'PidiQ 1 
'"'■/"•'• ( -^-Px&QidiQ'Pi - YPxd'QAQtpi + aPxQ^Q 1 ?! - ^P^&Q&Q 1 "ft 



dtd 3 xa 3 



^^d'QAQ 1 ^ - ^d'QjdiQ 1 ft + PxQid'diQ 1 ?! ) , 



(4.5) 



3 This reduces to, for example, eq. (53) of pm where a canonical kinetic term was considered, as expected. 



where we have used the fact that P^xf(HP^x) is order 0(e), and we have neglected the last term in 
the second line in fl4.5|) . Now (fO|) can be used to evaluate j3\ in the last line in (|4.5| ). Taking time 
derivative to both sides in ( f4.4| ) and keeping only the leading-order terms we have 



d 2 /?! = -a'P^iQ 1 



a 2 P 



x 



2H 



W ■ 



(4.6) 



As emphasized in [36] that the above equation cannot be inserted directly into the action ( |4.3| ) since 
it contains Q and therefore will change the order of the equations of motion for Q . However, it is 
very convenient to make use of the equations of motion derived from the second-order theory ( 3.12| ) to 
eliminate Q in ( |4.6| ). To the lowest-order of slow- varying parameters, the second-order action ( |3,12| ) 
reduces to the form 



5(2) 



X - I dtd 3 xa 3 P jX G IJ (q 7 Q j - d i Q I d l Q J + 0(e) 



and gives 



5L 



SQ 1 



a 3 P, x ( -3HQ! -Qi + d l diQi ) + 0(e) . 



(4.7) 



(4.8) 



This vanishes when the perturbations Q solve the free Gaussian theory, but (SL/SQ )|i will be non- 
zero when considering the full theory and of course the third-order interacting theory. Solving Q 
from (fljj) and inserting them into ( t4.6|) , it follows that 



a 2 /?! 



^iQ 1 - g^V 



l 



5L 



2Ha 5Q 1 



0(e) 



(4.9) 



Substituting (4.9) into (4.5) and performing a lot of integrals by parts, we finally get the equivalent 
third-order action 

-,3 / ! 



S^ 



did'.v { y I^PxxGu^k + ~P,xxx<j>i<>.i<>i< ! Q 7 Q J Q X 



1 



P 2 



\IAJr\K 



+ a 6 -G u Pxk - -tttGu^k Q j Q j Q 



a 3 P 



2 
xx 



AH 



(4.10) 



a 



W/^I : 



Gijcj) K d*Q 1 d i Q J Q" - —G I jP XK d l Q 1 diQ J Q 



iJr>K 



a 3 P 2 x 
2H 



In2„J*-2AK 



G IJ ^ K Q i d z Q J d^Q n + 



2 
8L 



SQ 1 



F\Q) 



with 



Fl (Q) = T§^ ( d- 2 (Qjd 2 Q J ) - ^QjQ J 



(4.11) 



1 
AH T \~ w ~ ■" ' 2 ( 

The last term t-^j F 1 \Q) in ( [4.10|) which is proportional to (<5L/<5Q 7 )|i can be absorbed by a 

fields redefinition of Q into new fields Q , as in the single- field case. It can be shown easily that the 
appropriate fields redefinition is 



(}' Q 1 - F 2 (Q) = Q 2 + ^ T QjQ J ~ ^¥d-\Qjd 2 Q' 



(4.12) 



Such a fields redefinition where F 1 are quadratic in Q 1 , has no effect on any of the 0(Q 3 ) terms in 
the third-order action ( 4. ICQ , and thus we may simply replace Q with Q 1 there. On the other hand, 
the fields redefinition indeed modifies the quadratic part of the action, i.e. the Gaussian action ( |3.12| ), 
which transforms as 



S^[Q]^S^[Q]- dtd 3 x 



8L 



8Q 1 



F\Q) , 



(4.13) 



the second term here cancels the last term in ( |4.10| ) exactly, which is proportional to the first-order 
equations of motion ^j 



5. Calculating the Three Point Function 



In this section, we proceed to calculate the scalar fields three-point functions (Q 1 (ki)Q J (k.2)Q K (IC3)) , 
with the third-order perturbative action ( f4.10| ) derived in the above section. The calculation of the 
three-point functions is standard, and thus we simply collect the final results here. 

1. Contribution from Q Q Q interaction. 

In conformal time rj, this interaction can be written as 



'I^fJ^tK 



with dimensionless coefficient 

f(l)IJK 



drjd x—f(i) IJK Q Q Q 



"FT f -1 

y f P,XxGlJ<j>K + ■^P,XXX<j>I<ftj4>K 



(5.1) 



2 



3/2 



G LJ + 



2X 



ai 



ZIJ 



u 
-e-k 

€ 



(5.2) 



where we have expressed fmuK i n terms of slow-varying parameters defined in Section 2.2. After a 
standard calculation we find the contribution from this term as 



^■2(,2J +.Kr, 



jt r 1 

z(2vr) 3 <5 3 (k 1 + k 2 + W'lf j* n . (2 fc3) }_„ dri^r/AkWe 



+ perms + c.c. 



Hi 



4fcfCi fvn rCn 



(5.3) 



(2^( kl + k 2 + k3)/ ( y p3 ' y S perms 



where K = k\ + fc 2 + ^3> an d an asterisk '*' denotes that the corresponding quantities are evaluated 
at horizon crossing k = aH. Here "permutation" means total 6 ways of simultaneously rearranging 
the indices /, J and K and momenta fei, k 2 and k% (i.e. the index T is always tied to k\, and so on). 
2. Contribution from Q Q Q interaction. 



UK 



(27r)^(k 1+ k 2 + k 3 )/ 2 ^ 



2Hi 



(2) * P 3 x*TliVk!) V * ^ 2 



rC-i /Co >vi *^o f?3 



+ 



+ perms 



(5.4) 



with 



P 2 



f(2)I.JK = -jG L jP,XK ~ ' 2 Gij((>k 



PjX 

2 



3/ ' ^3n. . 6iJ 2 „ 

G/j ( — e ^ + 2e i^ + -rr- e ^,x + e K 



2e 



P 



x 



(5.5) 



3. Contribution from d l Q I diQ J Q K interaction. 
(27r) 3 <5 3 (k 1+ k 2 + k 3 )./f* 



2Ht f 2 f h + h 2k t k 2 1 

3) * n,a(2fc 3 ) (kl ' k2)h v^~ + ^ + k ] ! prrms 



with 



fl". 



f(3)IJK = —Z-PxxGij&k = ~ I -i- ) -GljER • 



Px\ 3/2 n 



4. Contribution from d l Q diQ Q interaction. 

2tff 



(27r)^(k 1 + k 2 +k 3 )// 4 ^ 



UK 



(4) * ^,IL(2«f) 



ki • k? 



r . . Y.i>j k ikj kik 2 k 3 
-K -\ 77 1 — — I + perms 



K 



K 2 



(5.6) 



(5.7) 



(5.8) 



with 



f(A)IJK 



5. Contribution from Q d 2 Q d~ 2 Q interaction. 

2H$ 



3/2 _ /3n_ „ m 2 _ 

Gij [ ~ e K + 21k + — — ex,x 



2e 



P 



x 



(27r)^(k 1 + k 2 + k 3 )/ ( 



UK 



1.2J.2 pp' 



with 



P 2 



f(5)IJK = -TrF7 G ij4>K = -2 -£- G/je^ 



p 



x 



3/2 



2P" r V 2 

6. Contribution from the fields redefinition Q >-^ Q 1 — F (Q). 

p ,x 
W 



(5.9) 



(5.10) 



(5.11) 



Redefinition Q J t-> Q 1 + Jf (/> 7 QjQ J gives 



(27r) d <5 d (ki + k 2 + k3 



p ,x xi i 



m 



P 3 xMiWt)\ 2 



,x* 



3/2 



JnJK u 3 



e J G JA A;f + perms 



(5.12) 



and redefinition Q i— ► Q^ 

-(2^) 3 ,5 3 (k 1 + k 2 + k 



—r d~ 2 (Qjd 2 Q J ) gives 

fff /Px* x3/2 



^c,n*(2*f) v 2 



eG^kxk^ + perms 



(5.13) 



Bring all these contributions together, after some simplifications, we find an expression for the 
scalar fields perturbations three-point correlation functions 



\K 



(Q 1 (ki)Q J (k 2 )Q A (k 3 )> = (27r)^(ki + k 2 + k 3 



Hi 



{2P,x?/ 2 \\pk 2 



-A 



UK 



(5.14) 



with 
A UK 



G IJ e K - 
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G U e K 



4k 2 k 2 k 2 (1 k 1 + k 2 2k 1 k 2 



k 2 k 2 



() _- + 2 k\k 2 ^ + 2k 2 ) + ^ _ fc| 



K 



e 



+ G "( 3-^ + 4^ + 6?* 



K 2 
12P" 2 



P 



a 



rv-i /Co rVl HjQ Aj^ 



A" 



K 2 



+ (ki • k 2 ) 



^i^hkj k x k 2 k z 



K 



JJ 

+ —e K 



2X 



U \ 4fc/C-| Kicy /Co 



+ perms 



H 2 e 2 ej K 3 
and Aijk defined in (|2.32|) is related to the above by 



1 Px ~ 
■Auk = -^J-^-AuK ■ 



K 2 



(5.15) 



(5.16) 



Here the Hubble parameter H, effective speed of sound c s , slow- varying parameters e 1 etc. and A are 
all evaluated at the time of horizon-crossing k ~ aH. In general, with an arbitrary c s , the adiabatic 
mode exits the horizon at c s k = aH, as addressed by many authors [45, g^]. In this paper we assume 
that u = l/c 2 s -1 is of 0(e). 

The above result fl5.14|) -( [5.16[ ) reduces to that of |3g], where the three-point functions of multi- 
field models with canonical kinetic term P = X — V have been investigated. Our result shows 
the dependence of non-Gaussianity on these slow- varying parameters. As argued in J27]], where a 
very general single field model was considered, our results shows that the final non-Gaussianity is 



proportional to these small slow- varying parameters, except a parameter A. However, from the last 
line in ( [5,15 ) it is easy to see that the large non-Gaussianities would arise in models with —jrr 3* 1 
during inflation. Note that the standard choice of kinetic term corresponds to the case A = 0. Models 
with large u, i.e. c s <Cl may be another source of large non-Gaussianity, as argued in p7f |, and it will 
be of great interest, we would like to leave this for a future investigation. 



6. Conclusion 



In this paper, we have calculated the three-point function of curvature perturbation at the horizon- 
crossing, which arises from a general multiple- field inflation model with the action ( |2.1| ). We use 
6N formalism to relate the curvature perturbation £ and the multiple scalar fields fluctuations. The 
result is given by ( |2.39| ) and ( |5.14 )-( |5.16D , where the momentum-dependent amplitude A is given 
by ( 5.15| )-( |5.16 ). This result includes a large class of inflationary models, and can be viewed as a 
multiple- field generalization of previous results considering k-inflation Jlq], ghost condensation O] 



and DBI-inflation [45, 46, 47]. In the case of standard canonical kinetic term, our result reduces 
to the previous known results fl36| , |73| , Pq [, where the method of Lagrangian formalism or the fields 
equations recently developed is used. Our formalism would be helpful to analyze the dependence 
of non-Gaussianities on the structure of the inflation models, and also be useful to study the non- 
Gaussianities in multi-field inflationary models which will be constructed in the future. 

In the presence of many light scalar fields coupled to Einstein gravity with a non-canonical kinetic 
term, we define some slow-varying parameters in order to control the theory and to find the solutions 
perturbatively in these small parameters. These parameters can be seen as a generalization of the 
standard slow-roll parameters. However, in this paper we do not ascertain the conditions under which 
a particular P(X,<p • ) would admit such a slow-varying limit. If this limit breaks done, the theory 
would become rather complicated and the calculation would be less clear. Moreover, as addressed in 



27 1 and etc., in the presence of a non-canonical kinetic term, the smallness of these slow- varying 



parameters do not imply that the inflation itself is slow-rolling. 

In multi-field inflation models, this microphysically-originated non-Gaussianities produced at the 
horizon-crossing provide the initial conditions for the superhorizon evolution of the non-Gaussianities 
afterwards [35]. Even though the "initial" non-Gaussianity is small, after the superhorizon evolution, 
the final "primordial" non-Gaussianity would be significantly large, as in the curvaton scenarios [p3, 

n, m, i3, m, iiEi, Hi- 



In this work, we studied the non-Gaussianities which are generated due to the non-linear relations 
between curvature and inflaton perturbations. In general, primordial density perturbations can be 
generated not only from the inflaton field(s) during inflation. There are also other mechanisms which 
may generate perturbations and also significant non-Gaussianities. Among such possibilities, for 
example, the curvaton mechanism [I 



55, 56, 



59, 72 



|], has been proposed and some of 
their observational consequences including issues of the non-Gaussianities have also been investigated. 
As addressed in the introduction, in the multiple-field scenarios, curvature perturbation can gener- 
ally evolve after the horizon-exiting. Detectable non-Gaussianities can be produced when the curvature 
perturbation is generated from the entropy perturbations at the end of inflation p7l 32], or during 
reheating process J38|, \ 



7C[. Moreover, in string inspired inflationary models, other string effects, 



for example, cosmic string effect should also be considered [74 



In this work, we studied the non-linearity of primordial perturbations measured by three-point 
functions. Indeed, the deviation from a Gaussian distribution may also be induced by higher-order 
correlation functions. It is interesting to go forward to find the exact fourth-order action of curvature 
perturbation and investigate the non-Gaussianities from the four-point functions and trispectra |p2j , 



Finally, in this work we only considered the case of c s ~ 1, it would be interesting to extend the 
current result to the case of an arbitrary, especially a small c s , which would be another source of large 
non-Gaussianities. On the other hand, the studies of multiple DBI-inflation [jH| , |46|| show that in the 
limit of c s <C 1, the multi- field effects are suppressed by c s , and thus the multi- field models reduce to 
a effective single-field model. We would like to leave this for a future investigation. 
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A. General Structure of the Expansion of the Action 

In order to expand the action (3.2) to the third-order of Q 1 , we need the expansions for P(X, (f) 1 ) and 



EijE % i — E 2 . Firstly, we expand P and X as 

p _ p(0) + p(l) + p(2) + p(3) + . . . s 

x = x^+x^+x^+x^ + ... , (A-1) 

where P' n ' and X^ n > are 0(Q n ) pieces of P and X respectively. From fl3.4| )-(p!q), we have 

i P 3P^ X i 

X (2) = -QiQ 1 - fad'PidiQ 1 - -jj-^jQtQ 1 + -^L^^jQiQJ _ — Gljdi Q I d i Q J , 

P 3P^ P^ Jf 



(A.2) 



Expanding P to the third-order, we have 

p« = Pxi* 1 ' + p 7 g 7 , 

p(2) = p xX (2) + lp xx(X (D)2 + ip /jQ V + P^jrWQ 7 , 
p( 3 ) = P X X (3) + P XI I (1) I (2) + P X/ X (2) Q 7 

+ \P t xxx{xVf + IPxxiixMfQ 1 + IPxijXMQtQJ + ip^^Q^ . 

b 2 2 o 

(A.3) 

The first integral in the action (3.2) can be expanded as 
S A = I dtd 3 xVhNP 

= f dtd 3 xa 3 (l + «! + •••) (P (0) + P (1) + P (2) + P (3) + • • • ) (A.4) 

= f dtd 3 xa 3 (P(°) + P« + ai p(°) + P 2 + ai P« + P^ + ai p( 2 ) + • • • ) . 
Therefore the second and the third-order pieces of Sa are 

4 2) = I dtd 3 xa 3 (P( 2 ) + ai P«) , (A.5) 

4 3) = / dtd 3 xa 3 (P® + ai p( 2 )) , (A.6) 



respectively, where _p( n )'s are given by ([/ 

Now we consider the second integral in (|3.2| ), which we denote as Sb- First we may expand 



EijEV - E 2 as 



n = -{EijE Lj - E 2 ) = n<°> + n« + n( 2 ) + 



3H 2 + -=-^A + —rdidjPxdidjPx - rr^ft)^ + 



V2a 4 



2a 4 



(A.7) 



then 



5b = / dttf'x— n 



dtd 3 xa 3 



n(°) + (n« - ai n(°)) + (n( 2 ) - ai n« + « 2 n(°)) + (- ai n( 2 ) + a 2 n« - «fn(°)) + 

(A.8) 
The second and third-order pieces of the second integral in the coupled action (ET^) now read 



sg> = ('dtd 3 xa 3 (nM - ai u^ + a 2 n(°)) , 



5® = / dtd 3 xa 3 {-a 1 Ii {2) + a?n (1) - afu 



(A.9) 
(A.10) 



where IF n )'s are given by 

Finally, the second and third-order pieces of the effective action are 

s(2) = S ^ + sg } = J dtd 3 xa 3 (p^ + ai pd) + n ( 2 ) - a 1 n( 1 ) + « 2 n(°) ' 

cM 3 xa 3 (P^ + ai pW - ai n( 2 ) + a 2 uW - a 3 rt°A 



S (3) = s (3) + g 



(3) 
B 



(A.11) 

(A.12) 
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